Evidence for superconducting gap nodes in the zone-centered hole bands of KFe2As2 from magnetic 

penetration-depth measurements 
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Among the iron-based pnictide superconductors the material KFe2As2 is unusual in that its Fermi surface 
does not consist of quasi-nested electron and hole pockets. Here we report measurements of the temperature 
dependent London penetration depth of very clean crystals of this compound with residual resistivity ratio 
> 1200. We show that the superfluid density at low temperatures exhibits a strong linear-in-temperature de- 
pendence which implies that there are line nodes in the energy gap on the large zone-centered hole sheets. The 
results indicate that KFe2As2 is an unconventional superconductor with strong electron correlations. 
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I. INTRODUCTION 

The discovery of high transition temperature (Tc) iron- 
based superconductors (IBS) raises fundamental questions 
about origin of superconductivity.' The microscopic pairing 
interactions which give rise to superconductivity are inti- 
mately related to the structure of the superconducting energy 
gap which may be probed experimentally by studying the na- 
ture of the low-energy quasiparticle excitations. In particular, 
the presence of nodes in the energy gap signals an unconven- 
tional pairing mechanism, and the position of the nodes can 
be a strong guide as the exact form of the pairing interaction 
Vkk'. 

Based on analyses of spin-fluctuation mediated pair- 
ing models, several different gap structures have been 
proposed. The rich variety of possible pairing states has 
as its origin the unusual multiband electronic structure of 
the IBS. In most IBS there are disconnected quasi-two- 
dimensional hole and electron Fermi-surface sheets. The for- 
mer are centered on the F point in the Brillouin zone and 
the latter at the zone cornerji. Strong scattering between the 
electron and hole sheets, corresponding to a wavevector q ^ 
(tt, tt) leads to a nodeless gap with sign change between the 
hole and electron sheets (nodeless s± state). ^'^ However, if in 
addition to this there is strong low q scattering this can stabi- 
lize a state with nodes in the electron and/or hole bands with 
either s or c?-wave symmetry. •'"^ c-axis Fermi-surface disper- 
sion could also generate horizontal line nodes. ^'^'^ 

Thus far experimental studies have given evidence for two 
distinct types of nodal structure in IBS:' one is nodeless but 
may have significant anisotropy'^'i and the other has line 
nodes. '^ '^ The former fully gapped state seems consistent 
with the nodeless s± state, but it may also be explained by 
the conventional s (s++) state. ^" For the latter case, strong 
evidence of nodal lines has been obtained from the penetra- 
tion depth and thermal conductivity measurements in the com- 
pensated metals LaFePO (T^ = 6K)''' '^ and BaFe2(As,P)2 
(Tc = 30K)ji^ In both materials, recent studie a'^i^''^^ suggest 
that the nodes are most likely located on the electron bands 



near the zone corner of the Brillouin Zone, and that the hole 
bands centered at the F point remains fully gapped. 

To uncover how the superconducting gap structure is related 
to the microscopic pairing in IBS, studies of KFe2As2 may be 
particularly instructive. This superconductor has a relatively 
low Tc — 3 K and is the end member of the (Bai_j;Ka;)Fe2As2 
series. Unlike most other IBS where the volumes of the elec- 
tron and hole sheets are roughly equal, in KFe2As2 the vol- 
umes differ by one electron per unit celli2^ This causes a sub- 




FIG. 1. (Color online). The Fermi surface of KFe2As2 (labels in 
correspond to Table|ll(. (a) 2D cross-sectional representation of our 
band-structure calculations (note the pillow surface from band 1 does 
not appear in this cut), (b) Schematic cross-section with shapes of the 
various pockets determined by dHvA^^ and ARPES"*' experiments, 
(c) 3D view of the DFT calculated Fermi surface with bands energies 
shifted to best fit the dHvA frequencies.— 
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stantial change in the Fermi surface topology. As shown in 
Fig.[T] the ubiquitous X centered electron sheets are replaced 
by small quasi-two-dimensional hole-like tubes which do not 
nest at all with the F-centered hole sheets. If nesting does 
play an important role in the superconductivity of the high- 
Tc IBS then the nature of the superconducting state should be 
quite different in KFe2As2. It is also intriguing that the elec- 
tronic specific heat is quite large (7 ~ 93mJ/K^mol) com- 
pared with other IBS,^"* suggesting the importance of elec- 
tron correlations in KFe2As2. In this paper, we report mea- 
surements of the temperature dependent penetration depth of 
KFe2As2 which shows that this material has well formed line 
nodes located in the F-centered hole bands. 



II. EXPERIMENTAL METHODS 

The temperature dependence of the magnetic penetration 
depth A(r) was measured using a self -resonant tunnel-diode 
oscillator which was mounted in a dilution refrigerator.'^ 
The sample (approximate dimensions 0.2 x 0.2 x 0.02 mm"^) is 
mounted on a sapphire rod, the other end of which is glued to a 
copper block on which a RUO2 thermometer is mounted. The 
sample and rod are placed inside a solenoid which forms part 
of the resonant tank circuit which operates at ^ 14 MHz. The 
RF field within this solenoid is estimated to be < 10^^ T so 
that the sample is always in the Meissner state. DC fields are 
screened to a similar level using a mu-metal can. Changes in 
the resonant frequency are directly proportional to changes in 
the magnetic penetration depth as the temperature of the sam- 
ple is varied. The calibration factor is determined from the 
geometry of the sample, and the total perturbation to the res- 
onant frequency due to the sample, found by withdrawing the 
sample from the coil at low temperature .-22, The sapphire sam- 
ple holder has a very small paramagnetic background signal 
which varies ^ 1/{T + 9) which corresponds to a change of 
^ 0.4 nm in A of our samples between 100 mK and 200 mK. 
This is around ten times smaller than the signal from the sam- 
ple and was subtracted. The ac magnetic field is applied par- 
allel to the c axis so that the shielding currents flow in the 
a6-plane. To avoid degradation of the crystals due to reaction 
with moisture in the air, we cleaved the crystals on all six sides 
while they were coated in a thick layer of grease. The mea- 
surements were done just after the cleavage without exposure 
in air. The relatively sharp superconducting transitions found 
in the frequency shift of the oscillator as well as in the specific 
heat measured after the penetration depth measurements [in- 
set of Fig.m indicate that our procedure does not reduce the 
sample quality. 

The single crystals were grown by a self-flux method which 
will be described in detail elsewhere.— The temperature de- 
pendence of dc resistivity p{T) was measured by the standard 
4-probe method. Au contacts were evaporated after Ar plasma 
cleaning of the surface, which give contact resistance less than 
1 n. The voltage contacts have finite widths which gives un- 
certainty of the absolute value of p up to ^ 22%, but the tem- 
perature dependence is not affected by this uncertainty. 

In order to determine the bulk homogeneity of the super- 



conductivity in our samples, specific heat measurements were 
performed on the same single crystal sample as was used for 
the penetration depth measurements. Because of the small 
size of these samples (mass ^ 4 /ig), a modulated temperature 
method was used.^^ Briefly, the sample is glued to a 10 /im 
diameter chromel-constantan thermocouple and heated with 
modulated light from a room temperature LED via an optical 
fibre. This method has a high sensitivity but has poor absolute 
accuracy, so the values are quoted in arbitrary units. 



III. RESULTS AND DISCUSSION 

Dc resistivity p{T) measurements show that our crys- 
tals are extremely clean with the residual resistivity ratio 
RRR = /9(300K)/po of 1280 [Figs.^a) and (b)]. The 
low-temperature normal-state p(T) follows the Fermi-liquid 
dependence po + AT'^ with A = 0.030(7) ^ilcm/K^, and 
the residual resistivity po is estimated by extrapolation. We 
note our data is not consistent with the non-Fermi-liquid T^'^ 
behavior reported by Dong et al.,^^ but is consistent with 
another previous report.^'' The obtained relatively large A 
value follows the Kadowaki-Woods relation A = aKw7^r^ 
with flKW ~ 10^^/iricm(Kmol/mJ)'^, indicating that strongly 




FIG. 2. (Color online), (a) Temperature dependence of in-plane 
resistivity p(T) in a single crystal of KFe2As2. (b) The same data 
plotted against below 10 K. 



3 



10^ 



o 

:! 10 



10" 



10" 




Nb,Sn 




• /DSrjRti'O^ 



T/T 

0.10 ' 0.15 



0.25 



10 



10 



10 



JO 



10 



Y (mJ/molKO 



FIG. 3. (Color online). Coefficient A vs Sommerfeld constant 7 
(Kadowaki- Woods plot'") for various superconductors showing the 
Fermi-liquid AT^ dependence of p(T). Blue circles are for s-wave 
superconductors, while red squares are for unconventional supercon- 
ductors with nodes in the gap."*' The lines represent A = aKw7^ 



with flKw ~ 
(dashed line). 



10"^ (solid line) and 4 x 10"^^incm(Kmol/mJ)^ 



correlated electrons with large mass are responsible for the 
Fermi-liquid dependence. As demonstrated in Fig. [3] 
KFe2As2 is at the edge of heavy-fermion superconductors 
with nodes in the energy gap which are widely believed to 
have unconventional mechanisms of superconductivity!^ 

Figure|4]shows the low-temperature variation of the change 
in the penetration depth AA(T) = A(T) — A(0) in two sam- 
ples. In both samples a strong T-linear dependence is ob- 
served over a wide temperature range of 0.1 < T/Tc < 0.25. 
Such a strong T-linear dependence is distinctly different from 
the exponential dependence expected in the fully-gapped su- 
perconducting state and is instead consistent with gap with 
well-developed line nodes. A linear temperature dependence 
can only be explained by the presence of line nodes. This 
is in contrast to power law behaviors with exponents close to 
2 which could be consistent either with nodal behavior in the 
dirty limit or with strong impurity scattering in e.g. the intrin- 
sically fully gapped s± state. Evidence for line nodes in this 
compound has also been found in NMR, specific heat^ and 
thermal conductivity^^ measurements. 

Deviations from the T-linear behavior of A (T) are observed 
at the lowest temperatures which are likely due to a finite zero- 
energy density of states created by a small amount of impurity 
scattering. When impurity scattering is present in supercon- 
ductors with line nodes, the low-temperature AA(T) changes 
from T to T^, which can be described by the empirical for- 
mula AA(r) cx TV(r + T*)^ A fit to this formula [solid 
Unes in Fig.gl gives T* w 0.3 K for sample 1 and T* « 0.5 K 
for sample 2, which indicates relatively small levels of disor- 




FIG. 4. (Color online). Change in the penetration depth AA(r) 
at low temperatures in two samples. The data are shifted vertically 
for clarity. Inset shows the frequency shift A/ of the oscillator con- 
taining sample 1 , and relative change in the specific heat divided by 
temperature AC /T in the same sample. 
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FIG. 5. (Color online). Normalized superfluid density at low tem- 
peratures follows a T-linear dependence (dashed line). We used 
A(0) — 260 nm estimated from the Fermi-surface parameters. The 
solid lines are fits to the empirical formula involving the impurity 
scattering in superconductors with line nodes (see text). Inset com- 
pares the normalized superfluid density data for the two samples. To 
account for the factor of two difference in the slope of AA(r) at low 
temperatures, we used the doubled A(0) value for sample 2. 



der in these crystals. 

Although the temperature dependence of A is consistent be- 
tween samples the absolute values of dX/dT differ by a factor 
two. Although our calibration procedure linking the measured 
frequency shift to the change in A has proved to be highly ac- 
curate 5-10%) for polycrystalline elemental test samples 
and high-Tc cuprate superconductors such as YBa2Cu307_5 
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TABLE I. Contributions of each band to tlie normalized superfluid density A'^ (0) / A'^ (T) evaluated from DFT band structure calculations as 
well as the dHvA measurements.^^ me is the free electron mass. 



From DFT calculations From dHvA 



Sheet 


# holes 


DOS (eV) 


Up (eV) 




name 


# holes 


me 




1 (pillow at Z) 


0.002 


0.055 


0.16 


0.4 










2 (inner tube at F) 


0.258 


1.11 


1.466 


32 


a 


0.17 


6 


31 


3 (middle tube at F) 


0.342 


1.48 


1.3959 


29 


c 


0.26 


13 


23 


4 (outer tube at F) 


0.390 


1.52 


1.5098 


34 


p 


0.48 


18 


31 


5 (tubes near X) 


0.009 


1.328 


0.55 


5 


s 


0.09 


7 


15 


total 


1.00 


5.494 


2.577 


100 




1.0 




100 



(Refs.l25l andlssh. in some cases where there is large surface 
roughness of the cut edges AA may be overestimated. We ob- 
tain almost identical temperature dependence of the normal- 
ized superfluid density A^(0)/A^(T) if we use the doubled 
A(0) value for sample 2 [inset of Fig.|5|. This indicates that 
only the calibration factor has a factor of two differences be- 
tween the two samples and the whole temperature dependence 
is quite reproducible. A similar effect was found in LaFePO. 
The lower value found for sample 1 is likely to be more repre- 
sentative of the intrinsic value although in LaFePO the values 
of AA found by scanning SQUID spectrometry"" were around 
a factor two smaller than our lowest estimate, This has im- 
plications for our calculations of the superfluid density as will 
be discussed below. 

To evaluate the normalized superfluid density 
nsiT)/nsiO) = A2(0)/A2(r), we need the value of 
A(0), which we cannot measure directly in our exper- 
iment. Recent small-angle neutron scattering (SANS) 
measurements-'^ estimate A(T = 55 mK) w 200 nm and 
/iSR measurements-'^ give X{Tc/2) « 280 nm which also 
suggests A(0) ^ 200 nm. These values are close to the 
A(0) values calculated from the Fermi surface parameters 
[see TablelU. Regardless of the choice of A(0) value within 
the uncertainties, the obtained temperature dependence of 
ns{T) shows a T-linear behavior over a even wider range 
of T than A(r) itself [Figs.E] and i). This is expected 
because the 1 — a{T/Tc) dependence of ns{T)/ns{0) gives 
AA(T)/A(0) ja{T/Tc) + fa2(r/Tc)2 + . . . which leads 
to slightly concave (superlinear) temperature dependence for 
the penetration depth as we observed. 

In order to proceed with a more quantitative analysis of 
our results we have estimated the contribution of each of the 
Fermi surface sheets to the total superfluid density. As a first 
approach we have calculated the band structure of KFe2As2 
using density functional theory (DFT) using the WIEN2k 
package^^ and the experimental lattice constants and internal 
positions,-*^ 4 x 10^ fc-points (in the full Brillouin zone) were 
used for the calculations of plasma frequencies Wp, Fermi sur- 
face volumes and sheet specific density of states (DOS) which 
are reported in TableU The Fermi surface topology and band 
masses are very similar to those reported previously^'' [Fig. 
[Tl. The calculated total DOS and ujp coiTespond to the Som- 
merfeld constant 7 13.0mJ/K-^mol and A(0) = 76.6 nm, 
respectively. The experimentally observed 7 = 93 mJ/K-^mol 
(Ref.i24i) implies a renormalization of 7.2 in the total den- 



sity of states at the Fermi level. Assuming that the superfluid 
density is renormalized by the same factor leads to A(0) be- 
ing increased to 205 nm. We note that in Galilean-invariant 
(translation-invariant) systems like liquid ■^He the Fermi liq- 
uid coiTections to A cancel in the zero temperature limit4£ 
However, in crystalline solids this cancellation is thought not 
to occur'*' and indeed in heavy-Fermion systems a reduction 
in superfluid density consistent with the thermodynamic mass 
enhancement is observed4^ 

Although the general features of this Fermi surface cal- 
culation are confirmed by angle-resolved photoemission 
spectroscopy (ARPES)''^ and de Haas-van Alphen (dHvA) 
measurements,^-' the exact size of the various sheets and 
their warping are not. Also, dHvA measurements show that 
the mass renormalization effects vary between the different 
sheets. So as a second approach we estimate the contribution 
of the various sheets to the superfluid density directly from the 
dHvA measurements, assuming each sheet is a simple two di- 
mensional cylinder. As the largest j3 sheet [see Fig.[T| was not 
observed by dHvA, we estimate its volume from the total hole 
number constraint and its mass by using the measured specific 
heat. Note that for sheets with more than one extremal dHvA 
orbit we have taken the average. Also we have ignored any 
possible contribution from the small pillow sheet (band 1). 
The contribution of each sheet to the superfluid density were 
then calculated using A^'^ = ii^e^nH /m* — w^/c^, where 
uh and m* are the hole density and effective mass for each 
sheet. From this we obtain a total superfluid density which 
coiTesponds to A(0) = 260 nm. The fact that this procedure 
and the direct calculation from the band structure (including 
renormalization) give values of A(0) which compare favorable 
to the direct measurements by /iSR and SANS gives us con- 
fidence in the accuracy of the result. A key result from this 
analysis is that the contribution of the e band near the X point 
is small. It contributes only up to ^ 15% of the total super- 
fluid density. 

The temperature dependence of the normalized superfluid 
density ps — A-^(0)/A^(r) is plotted in Fig.|6] Here we have 
used values of A(0) which encompass the above estimates as 
well as a factor of 2 larger value in case our values of AA(r) 
are overestimated because of remaining surface roughness in 
sample 1 . These plots show that the superfluid density at T = 
Tc/3 does not exceed ^ 75% of the zero-temperature value 
which indicates that the nodes are located on at least one of 
the F-centered bands which have large contributions to the 
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FIG. 6. (Color online). Normalized superfluid density obtained by 
using representative values of A (0) . The lines are the fits to the nodal- 
gap model with parameters indicated in the figure. 



superfluid density. 

In a nodal superconductor a rapid decrease of superfluid 
density with increasing temperature is indicative of a small 
value of the gap slope near the nodes /iAq — l/rfA/d^jnodo 
(in a simple d-wave model"^"* for T ^ Tc.: Ps{T) ~ 1 — 
41n2//iAo). The presence of multiple Fermi surface sheets 
complicates the analysis in the present case, but it is rea- 
sonable to approximately model the gaps and Fermi veloc- 
ity values on the different sheets by globally averaged val- 
ues of and A. We then approximate the variation of 
the gap with in-plane Fermi surface angle (/> by A((/)) = 
niin(/iAo(/!), Aq).^"^ We note that for fi — 2 this produces 
a very similar form of ps{T/Tc) to the more usual d-wave 
form A{(j)) ~ Aq cos(20). An alternative way of modelling 
the gap would be to add higher harmonics to this lowest or- 
der d-wave form. However, this would introduce more fit- 
ting parameters if more than one extra harmonic was required. 
In the weak-coupling limit, /i is the only free parameter in 
this model as the temperature dependent gap can be calcu- 
lated self-consistently4i Here however, we allow for possible 
strong-coupling corrections to Aq and leave this as a free pa- 
rameter, fixing the temperature dependence of A to its weak- 
coupling d-wave form. 

As shown in Fig.|6] we find this model fits our data well 
for all assumed values of A(0). There is strong covariance 
between the two parameters /i and Aq because the low T slope 
is determined only by the product /iAq whereas the variable 
ratio /i/Ao influences the higher T behavior only weakly. We 



find the data can also be reasonably well fitted by the inclusion 
of a second isotropic gap accounting for ~30% of the total, 
but the large reduction of ps at low temperatures cannot be 
reproduced if we assume that only the e tubes have line nodes 
as they only contribute ^ 15% to the total superfluid density. 
This reaffirms our conclusion that nodes must be present on 
the r centered hole bands, although there is a possibility that 
some of the sheets could be fully gapped. 

As demonstrated in Fig. [3] in superconductors with strong 
electron correlations where the A and 7 values are strongly 
enhanced. Cooper pairs with finite angular momentum (p, d- 
wave, etc.) are favorable as these states with small probabil- 
ities in small pair distances reduce the Coulomb repulsion^^ 
Our results indicate that KFe2As2 with 7 and A values com- 
parable to some heavy-fermion superconductors has well- 
developed line nodes in the zone-centered large bands. This 
implies that the electron correlations play an important role 
in this compound and that an electronic (non-phononic) pair- 
ing mechanism is needed to overcome the Coulomb repulsion. 
Since the F-centered sheets in KFe2As2 are relatively large, it 
is likely that the intra-band spin-fluctuations are important to 
create the sign change inside these bands, which gives rise to 
line nodes. Recent theoretical calculations"*^ reveal that the 
spin susceptibility in heavily hole doped system has relatively 
weak momentum dependence, which suggests the relative im- 
portance of scattering vectors other than q ~ (tt, tt). 



IV. CONCLUSIONS 

In summary, our results indicating line nodes in the F- 
centered bands in KFe2As2 are consistent with d-wave or 
horizontal nodal state. This is different both from the node- 
less states found in (Ba,K)Fe2As2 (Refs.lol andlTTh and from 
the nodal s state with nodes in the electron bands^"^ which 
may be the case in BaFe2(As,P)2 i'^'^''^^ How this difference is 
linked to the changes in the Fermi surface will be an important 
clue towards a microscopic mechanism of superconductivity 
in iron-based superconductors. 
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